Full 0{aew) electroweak corrections to e+e HHZ * 

Zhang Ren-You^, Ma Wen-Gan^'^, Chen Hui^, Sun Yan-Bin^, and Hon Hong-Sheng^ 
1 CCAST (World Laboratory), P.O.Box 8730, Beijing 100080, P.R.China 
^ Department of Modern Physics, University of Science and Technology 
of China (USTC), Hefei, Anhui 230027, P.R.China 



Abstract 

We calculate the full 0{aew) electroweak corrections to the Higgs pair production 
process e^e~ — > HHZ at an electron-positron linear collider in the standard model, 
and analyze the dependence of the Born cross section and the corrected cross section on 
the Higgs boson mass mn and the cm. energy -y/s. To see the origin of some of the 
large corrections clearly, we calculate the QED and genuine weak corrections separately. 
The numerical results show that the corrections significantly suppress or enhance the Born 
cross section, depending on the values of mH and ^/s. For the cm. energy — 500 GeV, 
which is the most favorable colliding energy for HHZ production with intermediate Higgs 
boson mass, the relative correction decreases from —5.3% to —11.5% as tjih increases from 
100 to 150 GeV. For the range of the cm. energy where the cross section is relatively 
large, the genuine weak relative correction is small, less than 5%. 
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I Introduction 



One of the most important goals of present and future colliders is to study the electroweak 
symmetry breaking mechanism and the origin of the masses of massive gauge bosons and 
fermions. As we know, within the Higgs mechanism [J the electroweak gauge fields and 
fundamental matter fields (quarks and leptons) acquire their masses through the interaction 
with a scalar field (Higgs field) which has an nonzero vacuum expectation value (VEV). And 
the self-interaction of the Higgs field induces the spontaneous breaking of the electroweak 
S\J{2)l U(l)y symmetry down to the electromegnetic U(l)£;Af symmetry. 

The present precise experimental data have shown an excellent agreement with the predic- 
tions of the standard model(SM) except for the Higgs sector These data strongly constrain 
the couplings of the gauge boson to fermions i^zff ^Wff')^ ™^ gauge self-couplings, 
but say little about the couplings of the Higgs boson to fermions (^Hff) and gauge bosons 
i^HZZ and Xrww)- Recent LEP2 experiment suggests that the Higgs boson has the mass 
with a lower bound of 114.4 GeV and a upper bound of 193 GeV at the 95% confidence level 
[HI [1]. People believe that with the help of future high energy colliders, such as the CERN 
Large Hadron Collider (LHC) and Linear Colliders (LC's), precise tests of the Higgs sector 
become possible. In the past few years, many theoretical works have been contributed to 
studying the Higgs boson productions and the properties of Higgs couplings at future high 
energy colliders [3] 0. 

To reconstruct the Higgs potential and verify the Higgs mechanism experimentally, not 
only the Yukawa couplings and the couplings of the Higgs boson to gauge bosons should 
be measured, but also the Higgs self-couplings which include the trilinear coupling Xhhh 
and the quartic coupling Xhhhh- These Higgs self-couplings can be probed directly only by 
multi-Higgs boson production. Due to the fact that the cross sections for three Higgs boson 
production processes are much smaller than those for Higgs boson pair production [Tj [Hj, 
the quartic Higgs self-coupling remains elusive. Recently, the Higgs boson pair production 
processes have been widely considered, and the cross sections for these processes in the SM 
have been evaluated at linear colliders and hadron colliders. The phenomenology calculations 
show that it would be extremely difficult to measure the Higgs self-coupling Xhhh at the 
LHC jHj, and e^e~ linear colliders, where the study of the e^e~ — > HHZ and HHvv can be 
performed with good accuracy, represent a possibly unique opportunity for performing the 
study of the trilinear Higgs self-coupling [^. For the center of mass (cm.) energy ^/s from 
500 GeV up to about 1 TeV, the HHZ production with intermediate Higgs boson mass is 
the most promising process among the various Higgs double-production processes, since its 
cross section is relatively large and all the final states can be identified without large missing 
momentum. When the cm. energy ^/s exceeds 1 TeV, the e^e~ HHvu process becomes 
sizeable, and it is possible to measure the trilinear Higgs self-coupling Xhhh by using this 
process. Therefore, in the first stage of a LC {^/s < 1 TeV), e'^e~ HHZ is the most 
promising channel to measure the Higgs self-coupling Xhhh- 

Although the cross section for e~^e~ — > HHZ with intermediate Higgs boson mass is only 
about 0.1 ~ 0.2 fb for ^/s < 1 TeV, the measurement of the Higgs self-coupling Xhhh through 
the process at e'^e" colliders can be significantly improved. For example, C. Castanier et al, 
conclude that a precision of about 10% on the total cross section for e^e~ HHZ can be 
achieved leading to a relative error on Xhhh of 18% with the help of high integrated luminosity 
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J £ = 2 ah~^ after performing the detailed simulations of signal and background processes 
at the TESLA[Tni. Other simulations demonstrate that the Higgs self-coupling Xhhh can be 
extracted more accurately by using a discriminating variable, namely the invariant mass of 
the HH system, and one can expect the high sensitivity of the triple Higgs self-coupling with 
an accuracy to 8% and better in multi-TeV e'^e~ collisions 11^. Therefore, to determine the 
Higgs self-couplings and reconstruct the Higgs potential, the theoretical prediction of the cross 
sections for e^e~ HHZ at a LC within per-cent accuracy is necessary. For this purpose, 
we investigate the e^e^ HHZ process at a LC in detail and present the calculation of the 
cross section for e~^e~ — > HHZ with the full 0{aew) electroweak corrections in the SM in 
this paper. 

The paper is organized as follows: in section^ we present the calculations of the Born 
cross section for e~^e~ — > HHZ and the full 0{aew) electroweak corrections to this process. 
The numerical results and discussions are presented in section ITTll In the last section, a short 
summary is given. 



II Calculations 

In this paper we adopt the 't Hooft-Feynman gauge of the SM. At the tree level, there are 
six Feynman diagrams relevant to the process e~^e~ — > HHZ ( shown in Fig.l). In Fig.l only 
the second Feynman diagram (Fig. 1(b)) contains a trilinear Higgs self-coupling vertex. In the 
SM the Higgs potential can be expressed as 

V = I^H^ + !^H^ + "I^H^ (2.1) 

where v = {^\/2Gf) ■ The trilinear Higgs self-coupling constant, Xhhh, can be derived 
from this potential directly. By using the Higgs self-coupling constant Xhhh = Smjj/v and 
the relevant Feynman rules for gauge interactions, we can obtain the tree level amplitude 
A^tree and the cross section o"tree for e'^e~ — > HHZ. 

The 0{aew) (one- loop level) virtual corrections to e+e~ HHZ can be expressed as 

(27r)^ f 

O-virtual = O-tree (^virtual = TH^r-r^ / d^>3y'R.e (A^tree-A4*irtual) ; (2-2) 

2piVsJ ^ 

spin 

where pi is the cm. momentum of the incoming positron, d$3 is the three-body phase space 
element, and the bar over summation recalls averaging over initial spins -A^ virtual is the 
amplitude of the virtual Feynman diagrams, including self-energy, vertex, box, pentagon and 
counterterm diagrams. All the Feynman diagrams and the relevant amplitudes are created 
by FeynArts 3 |13j . and the Feynman amplitudes are subsequently reduced by FormCalc 
|14j . Due to the fact that the Yukawa coupling of Higgs/Goldstone to fermions is proportional 
to the fermion mass, we neglect the contributions of the Feynman diagrams which involve 
H — e — e, — e — e, — e — or G~ — — e vertex. 

As we know, the contributions of the virtual diagrams contain both ultraviolet (UV) and 
infrared (IR) divergences, which can be regularized by extending the dimensions of spinor 
and spacetime manifolds to = 4 — 2e and giving the photon a fictitious mass m^, 
respectively. In this paper, we adopt the complete on-mass-shell (COMS) renormalization 
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scheme |16) to fix all the renormalization constants. All the tensor coefficients of the one- 
loop integrals can be calculated by using the reduction formulae presented in Refs. |17j and 
[TH] . As we expect, the UV divergence contributed by the loop diagrams can be cancelled by 
that contributed by the counterterm diagrams exactly, while the IR divergence still exists. 
Therefore, the 0{a^^) virtual cross section (Jvirtuai is independent of the UV regularization 
parameter Cuv = 1/e — 7_E + log47r, but still a function of the IR regularization parameter 

In order to get an IR finite cross section for e^e^ HHZ up to the order of 0{a'^^), we 
should consider the 0{aew) corrections to e~^e~ HHZ due to real photon emission. By 
using the general phase-space-slicing algorithm ^| , the contributions to the photon emission 
process e+e~ HHZj are divided into a soft and a hard contribution, 

cereal = O"soft + 0"hard = ^tree ((^soft + Kurd) , (2.3) 

where the "soft" and "hard" refer to the energy of the radiated photon. The energy 
of the radiated photon in the cm. frame is considered soft and hard if < AE and 
E^ > AE, respectively. Both (Tgoft and fihard depend on the arbitrary soft cutoff AE/E^, 
where Eh = y/sjl is the electron beam energy in the cm. frame, but the real cross section 
(Treai is cutoff independent. In our calculations the soft cutoff AE/E^, is set to be very small, 
therefore, the terms of order AE/E^ can be neglected and the soft contribution can be 
evaluated by using the soft photon approximation analytically ^B] HO] 

dCsoft = — dcJtree— ^ / — TTT ( — — ) . (2-4) 

2^1" J\k^\<AE 2A:!^ \pi -k^ P2 ■ J 

Here k^ = \J\k^\'^ + m^. k^ = (fc^, k^) is the four momentum of the radiated photon, and 
pi and p2 are the four momenta of e"*" and e~, respectively. As shown in Eq. (|2.4)) . the soft 
contribution has a IR singularity at = 0. The IR divergence from the soft contribution 
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cancels exactly that from the virtual corrections. Therefore, Cvirtuai+soft) the sum of the 
0[a.%^) virtual and soft cross sections, is independent of the infinitesimal photon mass m^. 

The hard contribution is UV and IR finite. It can be computed numerically by using 
standard Monte Carlo technique. In this paper, our computation of the hard contribution 
<7hard is performed with the help of CompHEP |21j . which is a package for evaluation of tree 
level Feynman diagrams and integration over multi-particle phase space by using the adaptive 
multi-dimensional integration program Vegas j22j . 

Up to the order of ©(Og^), the corrected cross section for e+e" HHZ is just the sum of 
the O(ag^) Born cross section <Ttree, the 0{af^) virtual cross section fTvirtuai and the 0{a'^^) 
real cross section cireah 

Ctot = ftree + ^virtual + f'"real = Ctree (1 + ^^tot) , (2-5) 

where 6tot = (^virtual + ^soit + ^hard is the full 0{aew) electroweak relative correction. As 
we expect, the corrected cross section, crtot, is independent of Cuv and m-y, since it doesn't 
contains any UV or IR singularity. 

To discuss the origin of some of the large corrections, we need to calculate the photonic 
(QED) corrections and the genuine weak corrections separately. The QED corrections com- 
prise two parts: the QED virtual corrections f^rtuai '^hich contributed by the loop diagrams 
with virtual photon exchange in the loop and the corresponding parts of the counterterms, 
and the real corrections cTreai- Therefore the QED relative correction (5e can be expressed as 

= ^^rtual + '^soft + 4ard, (2.6) 

where 5^^.^^a^l 

= ^^rtuai/^tree, and the genuine weak relative correction 5w is defined as 

= <5tot - <5e- (2.7) 

III Numerical results 

For the numerical calculation we use the following SM input parameters jl2j . 

me = 0.510998902 MeV, = 105.658357 MeV, = 1.77699 GeV, 

rriu = 66 MeV, rUc = 1.2 GeV, mt = 174.3 GeV, 

nid = 66 MeV, = 150 MeV, rub = 4.3 GeV, 

mw = 80.423 GeV, mz = 91.1876 GeV, a'^iO) = 137.03599976. (3.1) 

By using the relevant SM parameters listed above, we obtain 

V = {ViGfY^''^ = mw sinOw / ^/^m:^ = 250.356 GeV. (3.2) 

Besides these SM input parameters, six more input parameters should be given in the nu- 
merical calculation, which are the cm. energy ^/s, the Higgs mass mn, the mass parameter 
of dimensional regularization Q, the UV regularization parameter Cuvi the IR regularization 
parameter and the soft cutoff /S.E / E^. Since the corrected cross section utot is indepen- 
dent of the mass parameter of dimensional regularization, we set Q = 1 GeV in the numerical 
calculation. 
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x/i [GeV] 


mn [GeV] 


I 


O'virtual+soft [fb] 
II 


RES (ERR) 


500 


115 


-0.13108006089013 


- 0.13108006072333 


-0.1311(2) 




150 


-0.05089427858862 


- 0.05089427861858 


-0.0509(1) 


800 


115 


-0.10810404820086 


- 0.10810404820051 


-0.1081(3) 




150 


-0.08643256423041 


- 0.08643256423012 


- 0.0864(2) 


1000 


115 


-0.09123978783411 


- 0.09123978977518 


-0.0912(3) 




150 


-0.08051191441546 


- 0.08051191441530 


- 0.0805(2) 


2000 


115 


-0.04961299068878 


- 0.04961299068876 


- 0.0496(2) 




150 


-0.04848771307528 


- 0.04848771307526 


- 0.0485(2) 



Table 1: The 0{aiJ cross section Cvirtuai+soft for e~^e HHZ process for various Higgs 
boson mass (115 GeV and 150 GeV) and cm. energy values (500 GeV, 800 GeV, 1000 GeV 
and 2000 GeV). 

In Tabled we present some numerical results of the 0{af^) cross section (Tvirtuai+soft for 
e+e^ — > HHZ, where the soft cutoff AE/Ei^ and the UV regularization parameter Cuv ai'e 
set to be AE/Ef, = 1Q~^ and Cuv = 0, respectively. The middle two columns, labelled with 
I and II, are corresponding to the cases of rriy = W~^^ GeV and niy = 1 GeV, respectively. 
Although the Monte Carlo statistical error of the cross section cTvirtuai+soft is the order of 
10~^, we reserve the output numbers in columns I and II with 14 digits. By comparing the 
two columns of output numbers, we find that the results are stable over 8 digits when varying 
the fictitious photon mass niy from 10~^° to 1 GeV. Therefore, we draw a conclusion that 
the O(ag^) cross section iTvirtuai+soft is independent of within the statistical error. The 
results with 4 (or 3) significant digits and the corresponding Monte Carlo integration errors 
are presented in the last column which is labelled with RES (ERR). 

Analogously, the UV finiteness of (Tvirtuai can also be checked numerically. We find that 
the numerical results of ^virtual are stable over 7 digits when varying the UV regularization 
parameter Cuv from to 10^ for various ^/s and mn- For simplicity, we do not present these 
numerical results in this section. 

In Fig. 2 we present the 0{aew) relative correction to e^e^ HHZ as a function of 
the soft cutoff AE/Eb, assuming mn = 115 GeV and ^/s = 500 GeV. As shown in this 
figure, both ^virtuai+soft (= "^virtual + ^soft) and (5hard depend on the soft cutoff AE/Eb, but 
the full 0{aew) electroweak relative correction 5tot is cutoff independent. To show the cutoff 
independence more clearly, we present cxtot, the corrected cross section for e~^e~ HHZ 
which includes the full 0{aeuj) electroweak corrections, with the statistical errors from the 
Monte Carlo integration in Fig. 3. As shown in Fig. 3, a clear plateau is reached for AE/Eb in 
the range 10~^ — 10~^ and the corrected cross section utot is obviously independent of AE/Eb- 

Until now, we have checked the m^, Cuv and AE/Eb independence of the full 0{aew) 
electroweak corrections to e~^e~ — > HHZ. In the following calculation, m-y, Cuv and AE/Eb 
are fixed to be 10~^ GeV, and 10~^, respectively. 

In Fig. 4 we present the Born cross section atree and the corrected cross section utot for 
e+e" HHZ as functions of the Higgs boson mass tuh for ^/s = 500, 800 and 2000 GeV, 
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Figure 2: The 0{aeu)) relative correction to e^e HHZ as a function of the soft cutoff 
AE/Eb. 
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Figure 3: The dependence of the corrected cross section for e+e HHZ on the soft cutoff 
AE/Eb. 
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Figure 5: The full 0{aew) electroweak relative correction to e+e~ HHZ as a function of 
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[GeV] 


ruH [GeV] 


O-trco [fb] 


o-tot [fb] 


Stot [%] 




115 


0.17493(2) 


0.1629(2) 


-6.9(1) 


500 


150 


0.071834(6) 


0.06357(6) 


-11.50(7) 




200 


0.49611(3) • 10-3 


0.3329(2) • 10-3 


-32.90(4) 




115 


0.17428(2) 


0.1740(3) 


-0.2(1) 


600 


150 


0.10840(1) 


0.1041(1) 


-4.0(1) 




200 


0.031802(3) 


0.02935(2) 


-7.71(7) 




115 


0.15868(3) 


0.1632(3) 


2.8(2) 


700 


150 


0.11665(2) 


0.1155(2) 


-1.0(1) 




200 


0.058846(7) 


0.05665(7) 


-3.7(1) 




115 


0.14156(3) 


0.1471(3) 


3.9(2) 


800 


150 


0.11363(2) 


0.1135(2) 


-0.1(2) 




200 


0.07246(1) 


0.0705(1) 


-2.7(1) 




115 


0.11293(2) 


0.1168(3) 


3.4(3) 


1000 


150 


0.09890(2) 


0.0983(3) 


-0.6(2) 




200 


0.07790(1) 


0.0753(2) 


-3.3(2) 




115 


0.07119(2) 


0.0704(3) 


-1.1(4) 


1500 


150 


0.06684(2) 


0.0634(2) 


-5.1(3) 




200 


0.06165(1) 


0.0569(2) 


-7.7(3) 




115 


0.05021(1) 


0.0473(2) 


-5.8(4) 


2000 


150 


0.04812(1) 


0.0435(2) 


-9.6(4) 




200 


0.04630(1) 


0.0408(2) 


-11.9(4) 



Table 2: The Born cross section atree, the corrected cross section cjtot and the full 0{aew) 
electroweak relative correction 5tot for various Higgs boson mass and cm. energy values. 

respectively. As shown in this figure, each solid curve has two spikes at the vicinities of 
rriH = 2m VI/ and mn = 2mz, which just reflect the resonance effects at rriH = 2m\Y and 
mn = 2mz, respectively. For ^/s = 2000 GeV, both cjtree and crtot are insensitive to niH, 
decrease very slowly as the increment of mn from 100 to 220 GeV. In contrast to the case of 

= 2000 GeV, the cross sections are sensitive to rriH when ^/s = 500 GeV. They decrease 
rapidly to zero as tuh increases to about 204 GeV. 

To describe the full 0{aew) electroweak corrections to the Born cross section for e'^e" — > 
HHZ quantitatively, we plot the full 0{aew) relative correction Jtotj defined as 5tot = (<7tot — 
o"tree)/ctree, as a function of m/f in Fig.5. For ^/s = 500 GeV, which is the most favorable 
cm. energy for e^e~ — > HHZ with intermediate Higgs boson mass, the relative correction 
is negative in the range of 100 GeV < mn < 204 GeV. It decreases from —5.3% to —11.50% 
as niH increases from 100 to 150 GeV. Since the cross section near threshold is very small, 
the large relative correction in this region is phenomenologically insignificant. For ^/s = 
2000 GeV, the relative correction is also negative in the range of 100 GeV < niH < 220 GeV. 
It decreases from —4.6% to —12.5% as the increment of mn from 100 to 220GeV. For 
^/s = 800 GeV, the relative correction is positive when mn < 150 GeV. It varies from 5.2% 
to —3.7% as rriH running from 100 to 220 GeV. The numerical results of iTtree, crtot and 5tot 
for some typical values of niH and ^/s are presented in Table |21 

In Fig. 6 we depict the Born cross section cTtree and the corrected cross section cJtot for 
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Figure 7: The full 0{aew) electroweak relative correction to e+e~ HHZ as a function of 
the cm. energy. 
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e~^e HHZ as functions of the cm. energy y/s for uih = 115, 150 and 200 GeV, respec- 
tively. The cm. energy varies in the range 500 GeV < ^/s < 2000 GeV, which is accessible at 
future hnear colliders, such as TESLA NLC HH, JLC ^ and CERN CLIC 2B]. From 
this figure we find that both the Born cross section atree and the corrected cross section cjtot 
increase firstly, reach their maximal values, and then decrease with the increment of y/s. For 
rriH = 115 GeV, atree and o"tot reach their maximal values of about 0.178 and 0.174 fb at 
-^/s ~ 550 GeV and ^/s ~ 600 GeV, respectively. For tjih = 150 GeV cJtree and atot reach the 
maximums of about 0.117 and 0.116 fb at ^/s ~ 700 GeV, and for rriH = 200 GeV they reach 
about 0.078 and 0.076 fb at \/s ~ 950 GeV, respectively. 

The dependence of the full 0{aew) electroweak relative correction to e^e~ — > HHZ on 
the cm. energy ^/s is displayed in Fig. 7. As shown in this figure, the full 0{aew) electroweak 
corrections suppress the Born cross section in the cm. energy range of 500 GeV < y/s < 
2000 GeV for niH = 150 and 200 GeV, while enhance the Born cross section in the cm. 
energy range of 610 GeV < ^/s < 1360 GeV for uih = 115 GeV. The relative corrections can 
reach about —6.9% and —11.5% at ^/s = 500 GeV for mn = 115 and 150 GeV respectively. 
For niH = 200 GeV, the relative correction is large and can exceed —10%. It ranges from 
-2.7% to -11.9% as varying in the range of 800 GeV < \/i < 2000 GeV. We can see 
that in some parameter space the electroweak relative corrections are only few percent and 
might be below the achievable experimental accuracy. 
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Figure 8: The dependence of the QED relative correction Se and the genuine weak relative 
correction 6w to e~^e~ — > HHZ on the Higgs boson mass. 

In Fig. 8 we present the dependence of the QED relative correction 6e and the genuine 
weak relative correction dy/ to e~^e~ — > HHZ on the Higgs boson mass separately. From this 
figure we find that at threshold the genuine weak relative correction is a striking contrast 
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to the QED relative correction, and the contribution of the 0{aew) electroweak correction is 

overwhelmingly dominantcd by the QED correction. For = 500 GeV, the genuine weak 
relative correction strongly depends on the Higgs boson mass. It increases from —1.20% to 
4.56% as niH increases from 100 to 150 GeV. For ^/s = 2000 GeV, the genuine weak relative 
correction is insensitive to the Higgs boson mass. It is about —15% ~ —17% for mn in the 
range of 100 GeV < tuh < 220 GeV. 

The ^/s dependence of the QED relative correction Se and genuine weak relative correction 
(5w to e+e" — > HHZ are displayed in Fig.9. Together with Fig.6 we can see from this figure 
that for the range of ^/s where the cross section is relatively large, the QED relative correction 
is not too large and the genuine weak relative correction is less than 5%. For > 1000 GeV, 
the Higgs mass dependence of the genuine weak relative correction is small. As -^s increases 
to 2000 GeV, the genuine weak relative correction can reach about —15% ~ —18% for niH in 
the range of 115 GeV < tuh < 200 GeV. 
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Figure 9: The dependence of the QED relative correction and the genuine weak relative 
correction to e~^e~ HHZ on the cm. energy. 



IV Summary 

In this paper we calculate the full 0{aew) electroweak corrections to e+e~ — HHZ at a 
LC in the SM, and analyze the dependence of the Born cross section, the corrected cross 
section including full 0{aew) electroweak corrections and the relative correction on the mn 
and To see the origin of some of the large corrections clearly, we calculate the QED and 
genuine weak corrections separately. From the numerical results we find that the full O 
electroweak corrections significantly suppress or enhance the Born cross section, depending 
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on the Higgs boson mass and the cm. energy of a LC. Both the Born cross section and the 
corrected cross section are insensitive to the Higgs boson mass in the range of 100 GeV < 
tuh < 220 GeV for ^/s = 2000 GeV, but strongly related to the Higgs boson mass in the 
range of 100 GeV < niH < 204 GeV for y/s = 500 GeV. With our chosen yfs — mu parameter 
space in this paper, the relative corrections are a few percent generally, and can exceed —10% 
when \fs ~ 500 GeV and mn ~ 150 GeV. Therefore, the full Oip-ew) electroweak corrections 
should be taken into account in the precise experiment analysis. We should also mention that 
in some parameter space, where the cross section is sizeable, the total relative corrections are 
only few percent and thus probably might be below the achievable experimental accuracy. 

Note Added: As we were amending this manuscript, we became aware of a similar paper 
by G. Belanger, et al..|27j. They presented a numerical comparison in Table 3 of R,ef.|27j. 
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Figure Captions 



Figure 1 The tree level Feynman diagrams for e~^e — > HHZ. 

Figure 2 The 0{acw) relative correction to e"*"e~ — HHZ as a function of the soft cutoff 
AE/Eb. 

Figure 3 The dependence of the corrected cross section for e+e~ HHZ on the soft 
cutoff AE/Eb. 

Figure 4 The Born cross section and the corrected cross section for e~^e~ — > HHZ as 
functions of the Higgs boson mass. 

Figure 5 The full 0{aew) electroweak relative correction to e+e~ HHZ as a function 
of the Higgs boson mass. 

Figure 6 The Born cross section and the corrected cross section for e~^e~ — > HHZ as 
functions of the cm. energy. 

Figure 7 The full 0{aew) electroweak relative correction to e+e~ HHZ as a function 
of the cm. energy. 

Figure 8 The dependence of the QED relative correction and the genuine weak relative 
correction to e^e^ HHZ on the Higgs boson mass. 

Figure 9 The dependence of the QED relative correction and the genuine weak relative 
correction to e'^e~ — > HHZ on the cm. energy. 

Table 1 The 0{af^) cross section (Jvirtuai+soft for e'^e~ HHZ process for various Higgs 
boson mass (115 GeV and 150 GeV) and cm. energy values (500 GeV, 800 GeV, 1000 GeV 
and 2000 GeV). 

Table 2 The Born cross section utree) the corrected cross section crtot and the full 0{aew) 
electroweak relative correction (5tot for various Higgs boson mass and cm. energy values. 
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